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Mean values and differential distributions of event-shape variables have been 
studied in neutral current deep inelastic scattering using an integrated luminos- 
ity of 82.2 pb _1 collected with the ZEUS detector at HERA. The kinematic range 



is 80 < Q 2 < 20 480GeV 2 and 0.0024 < x < 0.6, where Q 2 is the virtuality of 
the exchanged boson and x is the Bjorken variable. The data are compared 
with a model based on a combination of next-to-leading-order QCD calcula- 
tions with next-to-leading-logarithm corrections and the Dokshitzer- Webber non- 
perturbative power corrections. The power-correction method provides a reason- 
able description of the data for all event-shape variables studied. Nevertheless, 
the lack of consistency of the determination of a s and of the non-perturbative 
parameter of the model, ao, suggests the importance of higher-order processes 
that are not yet included in the model. 
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1 Introduction 

The hadronic final states formed in e + e~ annihilation and in deep inelastic scattering 
(DIS) can be characterised by a number of variables that describe the shape of the event. 
The event shapes presented in this paper are infrared- and collinear-safe observables and 
can be calculated using perturbative QCD (pQCD). In some cases, the prediction in next- 
to-leading order and next-to-leading-logarithm (NLO+NLL) approximation is available. 

Precision tests of the pQCD predictions using the experimentally measured event shapes 
require a good understanding of non-perturbative effects, namely the hadronisation pro- 
cess, which describes the transition from partons to the experimentally observed hadrons. 
These non-perturbative corrections decrease as a power of Q, the square root of the vir- 
tually of the exchanged boson, and are therefore called power corrections. They can be 
parametrised as X p /Q p , where the scale \ p and exponent p depend on the shape variable; 
the exponent p can be predicted by perturbation theory [1-4]. The success of this sim- 
ple model in fitting the data [5] has initiated many studies; previously non-perturbative 
effects could be estimated only through the use of Monte Carlo (MC) models. 

In the formulation of the power-correction model by Dokshitzer and Webber [2,6-8], 
the shape variables are given by the sum of the perturbative and non-perturbative parts 
which depend only on two constants: the strong coupling, a s and an effective low-energy 
coupling, ao, which is universal to all event shapes. This formulation allows the extraction 
of a s and Oq from a fit to the data. 

Studies of event shapes at HERA have been already reported by the HI and ZEUS 
collaborations [9,10]. This paper extends with increased statistics the previous ZEUS 
measurement of the mean event-shape variables to the analysis of differential distributions 
and of two new shape variables. Measurements were performed in the Breit frame [11] in 
the kinematic range 0.0024 < x < 0.6, 80 < Q 2 < 20 480GeV 2 and 0.04 < y < 0.9. Here 
x is the Bjorken variable and y = Q 2 /sx, where s is the centre-of-mass energy squared of 
the ep system. Inclusion of the differential distributions allows an improved test of the 
validity of the power correction method. 



2 Event-shape variables 

The event-shape variables studied in this analysis are thrust, T, jet broadening, B, the 
invariant jet mass, M 2 , the C-parameter, the variable y^ (defined below) and the mo- 
mentum out of the event plane, A'out- Thrust measures the longitudinal collimation of a 
given hadronic system, while broadening measures the complementary aspect. These two 



parameters are specified relative to a chosen axis, denoted by a unit vector n . Thus: 



T 



Hi\Pi 



B = Slf^,* 1 , (2) 

where pi is the momentum of the final-state particle i. 

When n is the direction of the virtual-photon, thrust and broadening are denoted by T 7 
and _B 7 , respectively. Alternatively, both quantities may be measured with respect to the 
thrust axis, defined as that direction along which the thrust is maximised by a suitable 
choice of n . In this case, the thrust and broadening are denoted by T? and Bt- 

In the Born approximation, the final state consists of a single quark, and T 7 and Tt are 
unity Consequently, the shape variables (1 — T 7 ) and (1 — Tt) are employed so that 
non-zero values at the parton level are a direct indicator of higher-order QCD effects. 

The normalised jet invariant mass is defined by 



|2 

where Ei is the energy of the final-state particle %. 
The C-parameter is given by 



M 2 = ( ^f^ ^ Pt] , (3) 



c = ^ j ^!-^Z . (4) 



3 Ei ? -|P?IIPj|sin 2 

where 9y is the angle between two final-state particles, i and j. 

The shape variables in Eqs. (HJ-(HJ) are summed over the particles in the current hemi- 
sphere of the Breit frame. To ensure infrared safety, it is necessary to exclude events in 
which the energy in the current hemisphere is less than a certain limit, Su m . The value 
£nm = 0.25 ■ Q was used [12]. The analysis is based on event shapes calculated in the 
P-scheme, i.e. with particles assumed to have zero mass after boosting to the Breit frame. 

In addition, two variables, yi and i^ouT, referred to as two-jet variables, are considered. 
The quantity y 2 is defined as the value of the jet resolution cut parameter, y CVLt , in the 
hr jet algorithm [13], at which the transition from (2+1) to (1+1) jets takes place in a 
given event; here the first number refers to the current jet(s) and the second to the proton 
remnant. 



The variable describing the momentum out of the event plane, i^ouT, has been suggested 
for study [14] in events with a configuration at least as complex as (2+1) jets. The event 
plane is defined by the proton momentum P in the Breit frame and the unit vector n 
which enters the definition of thrust major: 

T M = max ' _>. — , (5) 

Lib* I 

with the additional condition P ■ n =0. 
The variable -Kout is given by 

KovT = Y,W% (6) 

i 

where p° nt is the component of momentum p^ of the hadron i perpendicular to the event 
plane. For leading-order (LO) (2+1) configurations, since both jets lie in the event plane, 
only non-perturbative effects contribute to -Kout- At higher orders of a s , perturbative 
effects will also contribute to .Kout- 

In contrast to the definitions in Eqs. (JI])-(jIJ) the sums in Eqs. (JHJ) and (JBJ) run over all 
particles in the Breit frame. 



3 Experimental set-up 

The data used in this analysis were collected during the 1998-2000 running period, when 
HERA operated with protons of energy E p = 920 GeV and electrons or positrons 1 of 
energy E e = 27.5 GeV, and correspond to an integrated luminosity of 82.2 ±1.9 pb _1 . A 
detailed description of the ZEUS detector can be found elsewhere [15, 16]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are measured in the central tracking detector (CTD) [17], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 
polar- angle 2 region 15° < 6 < 164°. The transverse momentum resolution for full-length 
tracks can be parameterised as <j{j>t)Ipt = 0.0058^ © 0.0065 © 0.0014/pr, with px in 



1 In the following, the term "electron" denotes generically both the electron (e ) and the positron (e + ). 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 



GeV. The tracking system was used to measure the interaction vertex with a typical 
resolution along (transverse to) the beam direction of 0.4 (0.1) cm and also to cross-check 
the energy scale of the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [18] covers 99.7% of the total 
solid angle and consists of three parts: the forward (FCAL), the barrel (BCAL) and 
the rear (RCAL) calorimeters. Each part is subdivided transversely into towers and 
longitudinally into one electromagnetic section and either one (in RCAL) or two (in BCAL 
and FCAL) hadronic sections. The smallest subdivision of the calorimeter is called a 
cell. Under test-beam conditions, the CAL single-particle relative energy resolutions were 
a(E)/E = 0.18/V^ for electrons and a{E)/E = 0.35/VE for hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e'yp. The 
resulting small-angle energetic photons were measured by the luminosity monitor [19], a 
lead-scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 

4 Kinematics and event selection 

A three-level trigger system was used to select events online [16,20]. Neutral current DIS 
events were selected by requiring that a scattered electron candidate with an energy more 
than 4 GeV was measured in the CAL [21]. 

The offline kinematic variables Q 2 , x and y were reconstructed using the double angle 
(DA) method [22]. For offline selection the electron (e) and the Jacquet-Blondel (JB) [23] 
methods were also used. 

The offline selection of DIS events was based on the following requirements: 

• E' e >10 GeV, where E' e is the scattered electron energy after correction for energy loss 
in inactive material in front of the CAL, to achieve a high-purity sample of DIS events; 

• y e <0.9, where y e is y as reconstructed by the electron method, to reduce the photo- 
production background; 

• 2/jb>0.04, where t/jb is y reconstructed by the JB method, to ensure sufficient accuracy 
for the DA reconstruction of Q 2 ; 

• 38<£<60 GeV, where 5 = ^2 i (E — P z )i and the sum runs over all CAL energy deposits. 
The lower cut removed background from photoproduction and events with large initial- 
state QED radiation. The upper cut removed cosmic-ray background. For events 
with forward electrons with 0f ab < 1 radian, where the Gf ab is the polar angle in 
the laboratory system, the 5 cut was tightened to 44 < 5 < 60 GeV, to reduce the 
contributions from electromagnetic deposits outside the CTD that are likely to be 
neutral pions wrongly identified as electrons; 



• |Z vtx |<50 cm, where Z vtx is the Z position of the reconstructed primary vertex, to 
select events consistent with ep collisions. 

The kinematic range of the analysis is: 

80 < Q 2 < 20 480GeV 2 , 0.0024 < x < 0.6 and 0.04 < y < 0.9. 

For each event, the reconstruction of the shape variables and jets was performed using a 
combination of track and CAL information, excluding the cells and the track associated 
with the scattered electron. The selected tracks and CAL clusters were treated as massless 
Energy Flow Objects (EFOs) [24]. The minimum transverse momentum, p?, of each EFO 
was required to be greater than 0.15 GeV. 

The variables M 2 , C, T, and B were reconstructed only using objects in the current region 
of the Breit frame, with the following additional requirements: 

• number of EFOs (hadrons, in the case of theoretical calculations) in the current region 
of the Breit frame > 2; 

• Kab ! < 1-75, where 7/i a b is the pseudorapidity of an EFO as measured in the labora- 
tory frame. 

Jets were reconstructed using the kx cluster algorithm [13] in the longitudinally invariant 
inclusive mode [25]. The jet search was conducted in the entire Breit frame. For the 
%)2 variable, at least two EFOs (hadrons) had to be found in the Breit frame. Since the 
proton remnants were explicitly treated by the jet algorithm, all hadrons from the current 
and target hemispheres of the Breit frame were considered. 

The -ft'ouT variable was reconstructed in the entire Breit frame, with the following cuts 
required by theory [14]: ?7 Brcit a r < 3, to remove the proton remnants, and y 2 > 0.1, 
to avoid small values of Tm- In addition, l^b ! < 2.2 was required to select a region of 
well understood acceptance. 



5 Monte Carlo simulation 

A Monte Carlo event simulation was used to correct the data for acceptance and resolution 
effects. The detector simulation was performed with the Geant 3.13 program [26]. 

Neutral current DIS events were generated using the Djangoh 1.1 package [27], com- 
bining the Lepto 6.5.1 [28] generator with the Heracles 4.6.1 program [29], which 
incorporates first-order electroweak corrections. The parton cascade was modelled with 
the colour-dipole model (CDM), using the Ariadne 4.08 [30] program. In this model, 
coherence effects are implicitly included in the formalism of the parton cascade. The Lund 



string-fragmentation model [31], as implemented in Jetset 7.4 [32,33], was used for the 
hadronisation phase. 

Additional samples were generated with the Herwig 5.9 program [34], which does not 
apply electroweak radiative corrections. The coherence effects in the final-state cascade 
are included by angular ordering of successive parton emissions, and a cluster model is 
used for the hadronisation [35]. Events were also generated using the MEPS option of 
Lepto within Djangoh, which subsequently uses a parton showering model similar to 
Herwig. 

For Ariadne, the default parameters were used. The Lepto simulation was run with 
soft-colour interactions turned off, and Herwig was tuned 3 to give closer agreement 
with the measured shape variables at low Q; the CTEQ4D [36] parameterisations of the 
proton parton distribution functions (PDFs) were used. The MC event samples were 
passed through reconstruction and selection procedures identical to those of the data. 
The set of MCs used here ensures that the influence of both the parton level (Ariadne 
versus Lepto, Herwig) and the fragmentation (Herwig versus Ariadne, Lepto) on 
the systematic uncertainties is included. 

The generated distributions include the products of strong and electromagnetic decays, 
together with Kg and A decays, but exclude the decay products of weakly decaying 
particles with lifetime greater than 3 x 10~ 10 s. 

6 QCD calculations 

6.1 Perturbative QCD calculations 

The mean values and differential event-shape distributions were analysed using different 
perturbative QCD calculations. 

For the mean event shapes, NLO QCD calculations have been performed using the pro- 
grams DISASTER++ [37] and DISENT [38], which give parton-level distributions. To 
determine the theoretical a s dependence of the variables, both programs were run with 
the CTEQ4A proton PDFs with five a s sets [36] . The mean value of each shape variable 
was found to be linearly dependent on a s (Mz) in the range 0.110—0.122. The calcula- 
tions were performed with the renormalisation and factorisation scales [ir = xrQ and 
Up = XpQ, respectively, where for the central analysis xr and xp were set to 1. 

Infrared and collinear safety ensures that the mean values of event shapes can be computed 
with fixed order calculations [1,2]. However, in order to describe the differential distribu- 



Thc parameter PSPLT was set equal to 1.8; otherwise default parameters were used. 



tions in the phase space region where the perturbative radiation is suppressed (region of 
small values of the shape parameters), large logarithmic terms must be resummed. 

To obtain the theoretical predictions for the differential distributions, DISASTER+ + 
events were generated using the DISPATCH [39] program with the MRST99 [40] PDFs. 
The final predictions of the differential distributions, combining the NLO and NLL calcu- 
lations as well as the power corrections (see Section T6.2J) . were made using the DISRE- 
SUM [39] package described below. 

To calculate the perturbative part of the differential distribution, (da/dV)pT, where V is 
the event-shape variable, DISRESUM matches the NLL resummed perturbative calcu- 
lation of the differential distribution to the corresponding NLO distribution. The details 
of the resummed calculations depend on the type of shape variable i.e global, T 7 and 
B 7 , or non-global, M, C and TV [12,39,41-43]. Three possible types of matching were 
investigated: logR matching, similar to that used in e + e~ annihilation analyses, M and 
M2 matchings. The last two were specifically introduced for DIS processes [41]. In ad- 
dition, a modified matching technique can be used for the three types of matching. The 
modification to the matching ensures firstly that the integrated cross section has the cor- 
rect upper limit at V — Knax, where V max is the maximum of the distribution [41], and 
secondly that, if the fixed order distribution goes smoothly to zero at the upper limit, the 
matched-resummed distribution has similar behaviour. The modification requires that 
the ln(l/V) terms in the resummation are replaced by expressions of the form 



V 



1 



V \V r 



max 



(7) 



In addition, to ensure the correct upper limit to the distribution after non-perturbative 
corrections, as discussed in Section RT21 the shift to the distribution is multiplied by 

(8) 



The resummation can be expressed in terms of a rescaled variable, 1/xlV, instead of 1/V, 
where Xl is a logarithmic rescaling factor [43] . The values of p, p s and xl were set by 
default to 1, 2 and 1, respectively, and were varied, as explained in Section^ to estimate 
the theoretical uncertainties of the method. 



6.2 Non-perturbative QCD calculation: power corrections 

Before the data are compared to the pQCD predictions, the latter require correction for 
the effects of hadronisation. Dokshitzer and Webber calculated power corrections to the 



event-shape variables in e + e~ annihilation, assuming an infrared-regular behaviour of the 
effective coupling, a c s [2,6-8]. The technique was subsequently applied to the case of 
DIS [44] and has been used here. 

In this approach, a constant, ao, is introduced, which is independent of the choice of 
the shape variable. This constant is defined as the first moment of the effective strong 
coupling below the scale hi and is given by: 

ajj(fii) = — / a cS (n)dn, (9) 

N J 

o 

where \xi corresponds to the lower limit where the perturbative approach is valid. This is 
taken to be 2 GeV, as in the previous analyses [5,9, 10,45,46]. 

The theoretical prediction for the mean values of an event-shape variable, denoted by 
(V), is then given by 

(V) (a s ,M) = (V) PT (a.) + (V) pow (a s ,o^), (10) 

where (V) PT is calculated using the NLO QCD calculation, and (V) is the power 
correction. The power correction is given by 



4 MA l 



00 P ow = «v^7^- (ii) 



The values of ay for (1 — T T ), (1 — T 7 ), C and M 2 are respectively 2, 2, 3ir and 1. For 
By and B^ more complex expressions were used [8,41]. 

The variable M. is the 'Milan factor' of value 1.49 [47], which takes into account two-loop 
corrections; it has a relative uncertainty of about ±20%, due to three- and higher-loop 
effects. The term A\ is given by: 



A x = — Hi 



«o - <x a (fi R ) ~^:( lo s(— )+^- + 1 ) «« W 



(12) 



where Cp — |, K — ^ — ^ — 9^/5 Po — H ~ f-^/ an d -^/i the number of active flavours, 
is taken to be five. 

The power-corrected differential distributions are given by: 

>) = ^-W->> (13) 

where da/dV-p^ is calculated as described in the previous section. 



7 Analysis method 

The event shapes were evaluated for event samples in selected bins of x and Q 2 . The 
choice of the bin sizes [48] was motivated by the need to have good statistics and keeping 
the migrations, both between bins, and from the current to the target region within each 
bin, small. The kinematic bin boundaries are listed in Tabled 

The predictions that combine the pQCD calculations and the power corrections are fitted 
to the measured mean and differential distributions, with the exceptions of 2/2 and -Kqut, 
to extract the (ao, a s ) values. The theoretical predictions for 2/2 and Kout are not yet 
available, therefore, no attempt is made to extract (oo, a s ) from these variables, but they 
are compared to the NLO QCD calculations and MC predictions at parton and hadron 
levels. 

Separate x 2 ~^s to the mean values as a function of Q and to the differential distributions 
in bins of Q 2 are performed for each variable. The distributions, when calculated to NLO, 
diverge at small values of the shape variable. The divergence is removed when evaluating 
the integral to determine the mean values. Consequently for the mean values, the fixed- 
order NLO prediction is used combined with the power correction according to Eq. (fTU|). 
For the differential distributions, the divergence is removed by matching NLO to NLL 
using DISRESUM. The distributions are corrected for hadronisation following Eq. (JT5J) . 
For each observable, the fit was performed with a s (Mz) and a^ taken as free parameters. 

The fits to both the mean values and the differential distributions were made using the 
Hessian method [49] which uses a full error matrix that includes correlated off-diagonal 
terms due to the systematic uncertainties. Therefore, the statistical and systematic un- 
certainties are not quoted separately and appear as one 'fit error' in all the tables. 

The mean values of the event shapes were evaluated over the full measured kinematic 
range. The range used in the fits to the differential distribution has been defined in- 
dividually for each shape variable and each Q 2 range. The ranges are limited by the 
requirements that the pQCD predictions should be well defined within the bin used in the 
fit and that the range used should not extend above the LO upper limit for the variable. 
The first requirement was based on the ratio (NLL + NLO + power corrections) /(NLL 
+NLO); bins were omitted at low values of the shape variable, where the ratio showed a 
rapid fall, indicating that the power correction is not well defined in this region. Also the 
range 0.8 < V/V max < 1 was excluded from the fit for 1 — T 7 , B^ and M, where the LO 
upper limit is equal to V max , to avoid the region where theoretical predictions are sensitive 
to the details of the matching between NLO and NLL calculations (where the matching 
modification discussed in Section loTI has a large effect). 

The final ranges are summarised in Table El Since the theoretical predictions for dif- 



ferential distributions are reliable only at high values of Q 2 , the fit was restricted to 
Q 2 > 320 GeV 2 . 



8 Corrections 

In each (x, Q 2 ) bin, the Ariadne MC was used to correct for the event acceptance and 
the acceptance in each bin of each event-shape variable. The acceptance is defined as the 
ratio of the number of reconstructed and selected events to the number of generated events 
in a given bin. The acceptance generally exceeds 70% for all bins, except at extremes of 
the Q 2 range and at low y. 

Agreement was found between the uncorrected data and the predictions of Ariadne 
throughout the entire kinematic range of each event-shape variable (see Section 10.1), 
thus confirming its suitability for the purpose of correcting the data. The data were also 
compared with the Herwig predictions; here the agreement with data was satisfactory 
but slightly worse than when using Ariadne. The correction factors were evaluated as 
the ratios of the generated to the observed values in each (x, Q 2 ) bin. The correction 
procedure accounts for event migration between (x, Q 2 ) intervals, QED radiative effects, 
EFO-reconstruction efficiency and energy resolution, acceptances in p? and 8, and EFO 
migration between the current and target regions. These correction factors are all within 
15% of unity, and the majority lie within 10% for the mean values of the shape variables. 
The correction factors for the differential distributions are typically within 20% of unity. 



9 Systematic uncertainties 



A detailed study of the sources contributing to the systematic uncertainties of the measure- 
ments has been performed. The main sources contributing to the systematic uncertainties 
are listed below: 

• the data were corrected using a different hadronisation and parton-shower model, 
namely Herwig or Lepto, instead of Ariadne; 

• the cut y e was changed from 0.9 to 0.8; 

• the cut on ?/jb was increased from 0.04 to 0.05; 

• the cut on 5 was tightened from 38(44) < 5 < 60 GeV to 40(46) < 5 < 60 GeV; 
the harder cut was used to estimate any residual uncertainties in the photoproduction 
background; 
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• the measured energies of clusters in the calorimeter were varied by ±3%, ±1% and 
±2% for the FCAL, BCAL and RCAL, respectively, corresponding to the uncertainties 
of the associated energy scales; 

• the EFO cuts on rjiab and pr > 150 MeV were tightened to |r/i a b| < 1.5 and pr > 
200 MeV; the cuts were also removed. 

The largest systematic uncertainty arose from the choice of Herwig as the hadronisation 
model. The other significant systematic was due to the ?7i a b selection. The remaining 
systematics were smaller than or similar to the statistical uncertainties. 

To estimate the theoretical uncertainties for both the mean values and the differential 
distributions, the renormalisation scale was varied by a factor of two, and studies were 
made of the effects of changes to \xi and to the Milan factor. To give an indication of the 
uncertainties due to mass effects, the data were reanalysed using the i?-scheme. For the 
mean values, the CTEQ4 PDFs were replaced by the MRST99 set. For the differential 
distributions, the additional parameters p and p s , that ensure the correct behaviour of the 
matching and shift, were varied as shown in Tables [7| and |H1 The logarithmic rescaling 
factor, XL-, was changed to 1.5 [43] and the CTEQ5 PDF was used instead of MRST99. 



10 Results 

10.1 Mean values 

The mean values of the event-shape variables are compared with the Ariadne predictions 
in Fig. Ql In general, there is a good agreement between data and MC. However the MC 
tends to overestimate the shape variables at low Q 2 , in particular M 2 . The Ariadne pre- 
dictions at the parton level are also shown. The difference between the hadron and parton 
level demonstrates the contribution from the hadronisation process, as implemented in 
Ariadne. It should be noted that the parton level of Ariadne, defined by the parton 
shower model, does not have a rigorous meaning in pQCD [50] and should be taken as 
indicative only. The structure in the theoretical distributions results from the different 
x-ranges associated with the Q 2 bins, see Table [0 

The mean values of the event-shape variables (1 — Tt), Bt, M 2 , C, (1 — T 7 ) and B 1 as 
a function of Q were fitted, by varying a s and ao, to the sum of an NLO term obtained 
from DISASTER+ + , plus the power correction as given by Eq. (jllj) . The data and fit 
results are shown in Fig. El For all variables the theory fits the data well. For (1— T T ), the 
best fit results in a negative power correction, whereas theory predicts a power correction 
equal to that found for (1 — Ty). 
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The extracted values (ao, as) are shown in Fig. El and in Tables El and E] The contours on 
the plot represent one standard deviation errors, corresponding to about 30% confidence 
level (CL), as well as the 95% CL regions based on the fit errors as calculated using 
the Hessian method. The theoretical uncertainties are not shown but are given in the 
tables, since they result in a correlated shift to all fit results. The current world average, 
a s (M z ) = 0.1182±0.0027 [51], is also shown. 

The ao values are in good agreement with those previously published [10], but somewhat 
lower than those obtained by the HI and e + e~ experiments. The values of a s obtained 
from fits to (1 — Tr), -B 7 , C and M 2 are roughly consistent with each other, but some- 
what above the world average value. However, Tables El and 0] show that the theoretical 
uncertainties are substantial and strongly correlated between variables. Fits to B? and 
(1 — T 7 ) give values of a s which are inconsistent with the values obtained with the other 
variables, as already observed in the earlier ZEUS measurement. 

For a s , the dominant uncertainty is that due to variation of the renormalisation scale. 
For oq, the variation in the Milan factor gives the largest uncertainty except in the case of 
B T The PDF uncertainty was evaluated by replacing the CTEQ4 PDFs by the MRST99 
set. With the exception of (1 — T 7 ), the changes in the fitted («o, ens) are of the order of 
the Hessian fit error. For (1 — T 7 ), the power correction becomes positive and the fitted 
values of a s (ao) change to 0.1285(0.3541), values that are in closer agreement with the 
other variables. If the model were robust, the fitted values of a s would be independent of 
fij. However a dependence on /xj is clearly evident in the tables. In view of these results, 
no attempt to extract combined values of (ao, as) from the mean event shapes was made. 

10.2 Differential distributions 

The differential distributions of the event-shape variables for Q 2 > 320 GeV 2 are compared 
to the predictions of Ariadne in Figs. 0] and El For all variables, Ariadne describes 
the data well. The parton level of Ariadne is also shown. The difference between the 
hadron and parton levels can be taken as illustrative of the hadronisation correction. 

The differential distributions for (1 — T 7 ), B 7 , M 2 , C and (1— TV), for which the theoretical 
predictions are available, have been fitted with NLL + NLO + PC calculations as shown 
in Figs. El and The solid (dashed) bars show the bins that were used (unused) in the 
fit as described in Section [3 

None of the three matching techniques discussed in Section 6.1 is strongly preferred theo- 
retically. Although the modification terms should be used to ensure the correct behaviour 
of the cross section, all options included in DISRESUM have been used. The results 
of fits using six different matching options are shown in Fig. El and Tables El and El The 
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X 2 of the fits does not depend significantly on the form of matching used. The M2mod 
matching has been chosen for this analysis in view of the minimal dispersion of a s and 
oTo for this type of matching. Tables |S] and |H1 show that the M2mod and Mmod matching 
techniques give fitted a s and ao values that agree, in general, within the Hessian fit errors; 
an exception is B 1 , which shows a five standard deviation shift in ao when the match- 
ing is changed from M2mod to Mmod. The logRmod matching gives larger systematic 
changes in the fitted a s (o7o), of the order of two (one) standard deviation. In all cases, the 
unmodified matching schemes, which are theoretically disfavoured, result in larger shifts 
than the corresponding modified matching. It can be concluded that matching-scheme 
uncertainties are approximately twice the Hessian fit errors. 

The results of the fit to the differential distributions using the M2mod matching scheme 
are summarised in Tables [7| and |HJ The model gives a good description of the differential 
distributions for the global variables (1 — T 7 ) and B 1 over a substantial range of the shape 
variables; the x 2 /dof of the fit is close to unity. For the non-global variables, (1 — T T ), 
C, M 2 , the fit is less good, with the x 2 /dof lying in the range two to four. The fitted 
a s values are consistent with the world average. With the exception of C, the Txq values 
are consistent with those obtained from the mean values. Figure El shows that, for the 
global variables, the fitted values of a s and Oq are consistent with being independent of 
the Q range. The non-global variables show a larger sensitivity to the Q range, possibly 
reflecting the poorer \ 2 of the fits. 

Tables [7| and |H] also give the theoretical uncertainties in the fitted (oo, as) values. For 
a s , the dominant theoretical uncertainties result from the renormalisation scale and the 
logarithmic rescale factor. The power factors in the modification terms also give rise to 
significant uncertainties for all variables except (1 — T 7 ). In contrast to the results found 
for the mean values, changes in the Milan factor and /// have no significant influence on 
the fitted a s . In general, all systematic uncertainties, with exception of the check on 
PDFs, for cJo are large compared to the fit errors and comparable in size to those of the 
mean fits. 

An estimate of the influence of the fit range is given in the two final lines of Tables [7| and |HJ 
This estimate was obtained by changing the fit range by half a bin at low values of the 
shape variable, where the influence of the NLL terms is greatest. For a s , the effect of the 
change is a few percent for (1 — T 7 ) and B 1 . In contrast, the fit values for the non- global 
variables are significantly dependent on the fit range. A comparison with the differential 
distributions measured by HI shows reasonable agreement with this analysis. However, 
it should be noted that the two analyses differ in the kinematic range of the fits as well 
as in many details of the fits. The value of «o, given by HI, agrees with the interval, 
approximately 0.4—0.5, obtained in this analysis. 

As in the case of the mean values, the fitted values of (oq, as) for the differential dis- 
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tributions are inconsistent with one another, with the non-global variables, M 2 , C and 
(1 — Tt), yielding a lower a s than the global variables (1 — T 7 ) and _B 7 , irrespective of the 
matching scheme used. The uncertainties due to the fit range and theoretical parameters 
preclude a meaningful determination of the average values for a s and «o from the fits to 
the differential distribution. 

10.3 Measurement of y<i and i^ouT 

As discussed in Section 2, the analyses of the variable y 2 and Kqut were made in the full 
phase space of the Breit frame, including both the current and target regions. In contrast 
to the variables discussed previously, the correction to y 2 is expected to fall as l/Q 2 - 
Although the general form of the correction is known, the theoretical calculations are 
not yet available. Consequently, no fit has been made but the data have been compared 
to Ariadne and NLO predictions. The distribution of y 2 and the mean of y 2 as a 
function of Q are shown in Figs. 110b and [TUb . respectively, together with the Ariadne 
predictions at the hadron and parton levels. The figures show that Ariadne describes 
the y 2 distribution for Q 2 > 320 GeV 2 well, but overestimates the means at lower Q. 

In Fig. ITUh . the y 2 distributions are compared with the NLO distribution from DISENT 
calculated using a s (M z ) = 0.116. Except at the lowest y 2 value for high Q, the NLO 
predictions describe the data well. In Fig. HUH , the mean values of y 2 are plotted as a 
function of Q and compared with the NLO predictions. The agreement with the NLO 
predictions is good over the entire range of Q. 

The -ft'ouT variable measures the momentum out of the event plane defined by two jets 
and thus depends on a 2 at lowest order. The data are compared to Ariadne predictions 
at the parton and hadron level in Fig. ^2 For the differential distribution, at the parton 
level, Ariadne agrees well with the tail of the Kqvt distribution but peaks at a lower 
value than the data. The hadron-level prediction, on the other hand, describes the data 
well everywhere, indicating the importance of hadronisation corrections to this variable. 
The mean value of -Kqut/Q agrees well with the expectation from Ariadne at the hadron 
level. In contrast, the parton-level predictions lie below the data, with a difference that 
decreases with Q, again indicating the importance of hadronisation effects. 

11 Summary 

Measurements have been made of mean values and differential distributions of the event- 
shape variables thrust T, broadening B, normalised jet mass M 2 , C-parameter, y 2 and 
Kqut using the ZEUS detector at HERA. The variables T and B were determined relative 
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to both the virtual photon axis and the thrust axis. The events were analysed in the 
Breit frame for the kinematic range 0.0024 < x < 0.6, 80 < Q 2 < 20 480GeV 2 and 
0.04 < y < 0.90. The data are well described by the Ariadne Monte Carlo model. 

The Q dependence of the mean event shapes T, B, M 2 and C, have been fitted to 
NLO calculations from perturbative QCD using the DISASTER-I — |- program together 
with the Dokshitzer- Webber non-perturbative power corrections, with the strong coupling 
a s (M z ) and the effective non-perturbative coupling a® as free parameters. 

Consistent values of a s are obtained for the shape variables (1 — Tt), -B 7 , M 2 and C, with 
Oo values that agree to within ±10%. For Bt, the «o value agrees with other variables, 
whereas a s does not. The variable (1 — T 7 ) gives a s and «o values that are inconsis- 
tent with the other variables and, in contrast to the other variables, that are sensitive 
to the parton density used in DISASTERH — h For all variables, the renormalisation 
uncertainties, the dominant theoretical uncertainty, are three to ten times larger than the 
experimental uncertainties. Also the \xi parameter used in the power corrections gives 
large uncertainties. These may be indications for the need for higher orders in the power 
corrections. 

The program DISRESUM together with NLO calculations from DISPATCH has been 
used to fit the differential distributions for the event-shape variables (1 — TV), M 2 , C, 
(1 — T 7 ) and B 1 for Q 2 > 320 GeV 2 . A reasonable description is obtained for all variables. 
The modified matching schemes give fitted values of a s that are consistent with the world 
average. With the exception of C, the values of a® are consistent with those found from 
the mean values and lie within the range 0.4—0.5. 

Comparison between the a s and «o from the fits to different variables show, however, that 
the results are not consistent within the experimental uncertainties. The renormalisation 
uncertainties are still large. There is a considerable sensitivity to small changes in the 
kinematic range of fits, indicating problems in the theoretical description of the data. 
Also the choice of matching scheme produces variations of the order of the experimental 
uncertainties. 

The power corrections for the variables 2/2 and Aout are not yet available. For 2/2, the 
data are well described by NLO calculations. The variable -Kqut is well described by 
Ariadne predictions at the hadron level. A comparison of Kqvt with parton and hadron 
level predictions of Ariadne indicates the need for substantial hadronisation corrections. 

In summary, the power-correction method provides a reasonable description of the data 
for all event-shape variables studied. Nevertheless, the lack of consistency of the a s and 
«o determinations obtained in deep inelastic scattering, for the mean values in particular, 
suggests the importance of higher-order processes that are not yet included in the model. 
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Bin 


Q 2 (GeV 2 ) 


X 






1 


80 - 160 


0.0024 - 0.010 






2 


160 - 320 


0.0024 - 0.010 






3 


320 - 640 


0.01 - 0.05 






4 


640 - 1280 


0.01 - 0.05 






5 


1280 - 2560 


0.025 - 0.150 






6 


2560 - 5120 


0.05 - 0.25 






7 


5120 - 10240 


0.06 - 0.40 






8 


10240 - 20480 


0.10 - 0.60 




Table 1: 


The kinematic boundaries of the bins in 


x and Q 2 . 


Q 2 ( GeV 2 ) 


1-: 


T T 


M 2 


c 


1-T 7 


Bry 


320 - 640 


0.1- 


0.3 


0.05-0.2 


0.3-0.7 


0.1-0.8 


0.15-0.4 


640 - 1280 


0.05- 


0.3 


0.025-0.2 


0.2-0.7 


0.1-0.8 


0.15-0.4 


1280 - 2560 


0.05- 


0.3 


0.025-0.2 


0.2-0.7 


0.1-0.8 


0.1-0.4 


2560 - 5120 


0.05- 


0.3 


0.025-0.2 


0.1-0.7 


0.1-0.8 


0.1-0.4 


5120 - 10240 


0.05- 


0.3 


0.025-0.2 


0.1-0.7 


0.1-0.8 


0.05-0.4 


10240 - 20480 


0.05- 


0.3 


0.025-0.2 


0.1-0.7 


0.1-0.8 


0.05 - 0.4 



Table 2: Ranges used for fits to the differential distributions. 



20 



Variable 


l-T T 


Bt 


M 2 


c 


1-T 7 


Bry 


ol s {M z ) 


0.1252 


0.1149 


0.1231 


0.1263 


0.1456 


0.1231 


Fit error 


±0.0010 


±0.0008 


±0.0010 


±0.0006 


±0.0035 


±0.0022 


X 2 /dof 


0.4150 


0.4873 


1.4003 


0.4127 


0.9725 


2.6992 


correlation 


-0.5337 


-0.5719 


-0.5275 


-0.1133 


-0.9257 


0.7610 


Xr = 0.5 


-0.0070 


-0.0068 


-0.0077 


-0.0072 


-0.0095 


-0.0062 


Xr = 2.0 


+0.0085 


+0.0065 


+0.0091 


+0.0088 


+0.0104 


+0.0067 


M = 1.19 


+0.0026 


+0.0019 


+0.0024 


+0.0030 


+0.0034 


+0.0012 


M = 1.79 


-0.0023 


-0.0017 


-0.0021 


-0.0026 


-0.0029 


-0.0011 


Hi = 1 GeV 


+0.0056 


+0.0039 


+0.0052 


+0.0067 


+0.0075 


+0.0024 


M = 4 GeV 


-0.0061 


-0.0047 


-0.0057 


-0.0070 


-0.0077 


-0.0032 


E-scheme 


+0.0046 


+0.0031 


+0.0036 


+0.0033 


+0.0033 


+0.0013 


PDF 


-0.0010 


-0.0005 


-0.0015 


-0.0008 


-0.0172 


-0.0024 



Table 3: Results for a s (Mz) from the fit to the mean values of the shape variables. 
The fit error is the total experimental error including both statistical and experi- 
mental systematic errors. The correlation coefficients are those between the fitted 
values of a s [Mz) andao (see Table^. The theoretical uncertainties (see text) are 
also shown. 



Variable 


l-T T 


Dj' 


M 2 


C 


1-T 7 


-B 7 


c^ (2 GeV) 


0.4622 


0.4349 


0.4184 


0.4122 


0.2309 


0.4352 


Fit error 


±0.0047 


±0.0044 


±0.0074 


±0.0030 


±0.0167 


±0.0044 


Xr, = 0.5 


+0.0105 


+0.0316 


+0.0239 


+0.0094 


+0.0339 


+0.1625 


x R = 2 


-0.0036 


-0.0089 


-0.0111 


-0.0039 


+0.0063 


-0.1030 


M = 1.19 


+0.0360 


+0.0343 


+0.0258 


+0.0232 


-0.0507 


+0.0272 


M = 1.79 


-0.0280 


-0.0294 


-0.0210 


-0.0198 


+0.0252 


-0.0201 


E-scheme 


+0.0157 


+0.0079 


+0.0120 


+0.0130 


-0.0043 


+0.0072 


PDF 


+0.0139 


+0.0113 


+0.0169 


+0.0129 


+0.1232 


+0.0032 



Table 4: Results for ttQ from the fit to the mean values of the shape variables. The 
fit error is the total experimental error including both statistical and experimental 
systematic errors. The theoretical uncertainties (see text) are also shown. 
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Variable 


1-T T 


M 2 


C 


1-T 7 


Bj 


a s (M z ) 


0.1151 


0.1158 


0.1176 


0.1227 


0.1226 


Fit error 


±0.0016 


±0.0013 


±0.0016 


±0.0012 


±0.0013 


Mmod 


-0.0009 


-0.0001 


±0.0007 


-0.0011 


-0.0020 


M 


-0.0020 


-0.0040 


-0.0033 


-0.0012 


-0.0031 


Ml 


-0.0038 


-0.0048 


-0.0011 


-0.0025 


-0.0036 


logRmod 


-0.0031 


-0.0032 


-0.0029 


-0.0029 


±0.0012 


logR 


-0.0045 


-0.0056 


-0.0026 


-0.0040 


-0.0054 



Table 5: Results for a s (Mz) from the fit to the differential distributions of the 
shape variables using the M2mod matching scheme. The fit error is the total ex- 
perimental error including both statistical and experimental systematic errors. The 
theoretical uncertainties due to the use of different matching schemes (see text) are 
shown. 



Variable 


l-T T 


M 2 


C 


1-T 7 


-B-y 


c^ (2 GeV) 


0.4173 


0.4650 


0.3358 


0.4820 


0.4268 


Fit error 


±0.0134 


±0.0100 


±0.0138 


±0.0138 


±0.0217 


Mmod 


-0.0012 


±0.0013 


±0.0042 


±0.0162 


±0.1048 


M 


±0.0041 


±0.0264 


±0.0496 


±0.0129 


±0.1329 


Ml 


±0.0087 


±0.0287 


±0.0657 


±0.0272 


±0.1506 


logRmod 


±0.0018 


±0.0130 


±0.0165 


±0.0159 


±0.0286 


logR 


±0.0107 


±0.0335 


±0.0650 


±0.0360 


±0.1707 



Table 6: Results for a® from the fit to the differential distributions of the shape 
variables using the M2mod matching scheme. The fit error is the total experimental 
error including both statistical and experimental systematic errors. The theoretical 
uncertainties due to the use of different matching schemes (see text) are shown. 
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Variable 


\-T T 


M 2 


c 


1-T 7 


-B-y 


a s (M z ) 


0.1151 


0.1158 


0.1176 


0.1227 


0.1226 


Fit error 


±0.0016 


±0.0013 


±0.0016 


±0.0012 


±0.0013 


X 2 /dof 


3.33 


2.46 


3.97 


0.74 


0.50 


correlation 


-0.72 


-0.72 


-0.64 


-0.69 


-0.75 


x R = 0.5 


-0.0023 


-0.0039 


-0.0039 


-0.0040 


-0.0028 


x R = 2 


+0.0054 


+0.0057 


+0.0051 


+0.0060 


+0.0049 


M = 1.19 


+0.0000 


-0.0000 


+0.0003 


+0.0001 


+0.0000 


M = 1.79 


-0.0000 


+0.0000 


-0.0002 


-0.0001 


-0.0001 


IH = 1 GeV 


+0.0001 


-0.0001 


+0.0004 


+0.0002 


+0.0001 


m = A GeV 


-0.0001 


+0.0000 


-0.0007 


-0.0004 


-0.0002 


x L = 1.5 


+0.0045 


+0.0046 


+0.0049 


+0.0024 


+0.0038 


x L = 0.67 


-0.0044 


-0.0048 


-0.0044 


-0.0034 


-0.0062 


p= 2.0 


-0.0029 


-0.0044 


-0.0069 


-0.0011 


-0.0038 


Ps = 1.0 


+0.0029 


+0.0031 


+0.0025 


+0.0016 


+0.0013 


E-scheme 


-0.0049 


+0.0033 


-0.0114 


+0.0009 


-0.0004 


PDF 


+0.0000 


+0.0003 


+0.0004 


+0.0009 


+0.0008 


-0.5 Bins 


-0.0143 


-0.0112 


-0.0066 


-0.0019 


-0.0022 


+0.5 Bins 


+0.0103 


+0.0073 


+0.0086 


+0.0039 


+0.0037 



Table 7: Results for a s (Mz) from the fit to the differential distributions of the 
shape variables. The fit error is the total experimental error including both sta- 
tistical and experimental systematic errors. The correlation coefficients are those 
between the fitted values of a s (M z ) andlx® (see Table\$). The theoretical uncer- 
tainties (see text) are also shown. 
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Variable 


l-T T 


M 2 


C 


1-T 7 


-Sy 


c^ (2 GeV) 


0.4173 


0.4650 


0.3358 


0.4820 


0.4268 


Fit error 


±0.0134 


±0.0100 


±0.0138 


±0.0138 


±0.0217 


x R = 0.5 


-0.0419 


-0.0344 


-0.0366 


-0.0444 


-0.0637 


x R = 2 


+0.0114 


+0.0215 


+0.0233 


+0.0284 


+0.0503 


M = 1.19 


+0.0453 


+0.0571 


+0.0211 


+0.0505 


+0.0394 


M = 1.79 


-0.0301 


-0.0377 


-0.0150 


-0.0349 


-0.0258 


x L = 1.5 


+0.0231 


+0.0308 


+0.0289 


+0.0442 


+0.0665 


x L = 0.67 


-0.0203 


-0.0230 


-0.0262 


-0.0385 


-0.0299 


p = 2.0 


+0.0042 


+0.0223 


+0.0424 


+0.0109 


+0.1010 


Ps = 1.0 


+0.0009 


+0.0029 


-0.0098 


-0.0014 


+0.0145 


E-scheme 


+0.0006 


+0.0026 


+0.0410 


+0.0107 


-0.0030 


PDF 


+0.0009 


-0.0019 


-0.0023 


-0.0035 


-0.0027 


-0.5 Bins 


+0.1473 


+0.1079 


+0.0751 


+0.0437 


+0.0514 


+0.5 Bins 


-0.0021 


-0.0203 


-0.1039 


+0.0348 


-0.0688 



Table 8: Results for a^ from the fit to the differential distributions of the shape 
variables. The fit error is the total experimental error including both statistical and 
experimental systematic errors. The theoretical uncertainties (see text) are also 
shown. 
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Figure 1: The mean values of event-shape variables as a function of Q. The 
solid and dashed lines are the hadron and parton level predictions from ARIADNE, 
respectively. 
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Figure 2: The mean values of event-shape variables as a function of Q. The 
solid lines are the results of the fit to the data of the predictions of the sum of NLO 
pQCD calculations from DISASTER+ + and the power corrections. The dashed 
lines are the DISASTERH — |- contribution to the fit alone. 
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Figure 3: Extracted parameter values for (a s , ao) from fits to the mean values 
of the shape variables. The vertical line and shaded area indicate the world average 
value of a s (M z ) [51]. 
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Figure 4: Differential distributions for the event shapes M 2 , C and 1 — Tt- 
The distributions are normalised such that n refers to the number of events in the 
(x, Q 2 ) bin after the Su m cut and N to the total number of events in the (x, Q 2 ) 
bin before the S\ im cut. The differential cross section has been scaled for clarity by 
factors 10 n ; where n= 12, 10, 8, 6, 4, 2 for < Q >= 21, 29, 42, 59, 82 and 113 
GeV, respectively. Predictions of ARIADNE at the hadron (solid lines) and parton 
(dashed lines) levels are shown. 
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Figure 5: Differential distributions for the event shapes 1 — T 7 and _B 7 . Other 
details as in Fig. ^| 
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Figure 6: Differential distributions for the event shapes M 2 , C and 1 — Tt- 
The distributions are normalised such that n refers to the number of events in the 
(x, Q 2 ) bin after the Su m cut and N to the total number of events in the (x, Q 2 ) 
bin before the S\ im cut. The differential cross section has been scaled for clarity by 
factors 10", where n= 12, 10, 8, 6, 4, 2 for < Q >= 21, 29, 42, 59, 82 and 113 
GeV, respectively. The solid(dashed) curves show the points used (omitted) in the 
fit to NLL resummed calculation matched to NLO plus power corrections. 
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Figure 7: Differential distributions for the event shapes 1 — T 7 and £? 7 . Other 
details as in Fig. E| 
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Figure 8: Extracted parameter values for (a s , oq) from fits to differential distri- 
bution of shape variables for the different matching techniques, which are discussed 
in SectionWIl 
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Figure 9: Extracted parameter values for (a s , ao) from fits to differential distri- 
bution of shape variables using M2mod matching for two different Q 2 intervals. 
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Figure 10: (a), (c) Differential distribution for y<x, scaled for clarity by different 
factors for each < Q > -value, and (b), (d) mean distribution of y2 versus Q com- 
pared with predictions from (a), (b) ARIADNE and (c), (d) NLO QCD calculations 
using DISENT. The ARIADNE predictions are shown for the hadron (solid) and 
parton (dashed) levels. The NLO calculation (shaded band) uses the CTEQ4A3 
proton PDFs. 
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Figure 11: (a) Differential distribution for Kovt/Q, scaled for clarity by dif- 
ferent factors for each < Q > value, and (b) mean distribution of -Kout/Q versus 
Q, compared with predictions from ARIADNE for the hadron (solid) and parton 
(dashed) levels. 
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